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Abstract 

We establish asymptotic formulas for the determinants of N x N 
Toeplitz + Hankel matrices Tn{4>) + Hn(<1>) as N goes to infinity for sin- 
gular generating functions cj> defined on the unit circle in the special case 
where <j> is even, i.e., where the Toeplitz + Hankel matrices are symmetric. 



1 Introduction 

In the theory of random matrices, for certain ensembles, one is led to consider the 
asymptotics of Fredholm operators of the form I + W + H where W is a finite 
and symmetric Wiener-Hopf operator and H is a finite Hankel operator fll3|| . 
This problem arises when investigating the probability distribution function 
of a random variable thought of as a function of the eigenvalues of a positive 
Hermitian random matrix. For general information about random matrix theory 
we refer the reader to jl3| and also to [j], 0, [t]] for more specific tie-ins to the 
random variable problem. 

The focus of this paper is to study the discrete analogue of this problem. 
This is not precisely the desired situation for those interested in random matrix 
theory. However, it is a natural starting place for cases where the random 
variable is discontinuous, since then the discrete nature of the computations 
make things a bit more accessible and the mathematical questions that arise 
are quite interesting in themselves. 

* ebasoriacalpoly.edu. Supported in part by NSF Grant DMS-9970879. 
ttehrhard@mathematik.tu-chemnitz.de. 
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The discrete analogue of this problem is to hnd an asymptotic expansion of 
the determinants of Toeplitz + Hankcl matrices 

M N (<p)=T N (cl>)+H N (cl)) (1) 

in the case where these matrices are symmetric. Here the N x N Toeplitz and 
Hankel matrices arc defined as usual by 

Tn{4>) = {(f>j-k)*^ , H N {(f>) = Oj+fc+O^Jo- (2) 

The entries 4> n are the Fourier coefficients 

<f> n = ±- <j>{e i9 )e- ine d9 (3) 
2tt Jo 

of a function cj> G L 1 (T) defined on the unit circle T. 

The matrices Mjv(0) are symmetric if and only if the function 4> is even, 
i.e., if (^(e 10 ) — <ft(e~ ie ). From the point of view of random matrix theory, one 
is particularly interested in the asymptotics of det Mn{4>) as iV — > oo in the 
case of even, piecewise continuous functions <p. This is related to the problem of 
finding the distribution function of a random variable that counts the number 
of eigenvalues of a random matrix that lie in an interval and to finding the 
distribution function for other random variables. See j| to see the connections 
between these problems. 

The problem of determining the asymptotics of the determinants of (not 
necessarily symmetric) matrices Mjv(0) has been studied intensively in a pre- 
vious paper [Q. For example, it was shown there that if <fi is continuous and 
sufficiently smooth, then the asymptotics are very similar to the ones given in 
the Strong Szego-Widom Limit Theorem. Indeed, it is only in the constant, or 
third order term that the answers differ. This is no surprise since if <fi is con- 
tinuous, then the Toeplitz operator is perturbed by a compact Hankel operator 
only. 

If, however, the symbol <f> is singular, then the problem is much harder 
to solve. In the case of Toeplitz determinants the answer is provided by the 
Fisher-Hartwig conjecture, which has been proved under certain smoothness as- 
sumptions in all the cases where it is expected to hold || . In Q an asymptotic 
formula for the determinants detM/v(0) was obtained for piecewise continuous 
functions (j>, but under the additional assumption that the function <fi does not 
possess a discontinuity at both a point on the unit circle and its complex con- 
jugate. In this case, the asymptotic formula shows that the asymptotics differ 
from the asymptotics of Toeplitz determinants only in the third order term, i.e., 
in the constant term, while the second order term is the same. 

Unfortunately, the additional assumption on the location of the discontinu- 
ities imposed in Q excludes all even, piecewise continuous functions. Hence the 
paper 0] does not answer the discrete analogue of the problem motivated by 
random matrix theory. 
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It is the purpose of this paper to solve this problem by establishing an 
asymptotic formula for determinants of matrices (|l|) for even piecewise contin- 
uous functions <fi. For such functions, it turns out that the asymptotics differs 
also in the second order term in comparison with the asymptotics of Toeplitz 
determinants. 

The paper is organized as follows. In Section ^ we will recall some of the 
results established in jl that are of relevance for this paper. In Section |^ 
we establish an identity which is the key for computing the asymptotics of 
det Mat (0) for <j> even. This identity can be formulated as follows: 



In this identity a can no longer be considered as a function, but has to be 
understood as a distribution, which satisfies certain properties. Moreover, a is 
here a certain concrete distribution. This identity was established first in ]12| 
in a formulation that is not based on distributions. The goal of Section || is to 
provide the necessary tools needed for dealing with distributions, in particular, 
to define a product between a and a in an appropriate way. Having done this, 
we are able to derive the distributional formulation of this identity from the 
original one. 

The identity (^) reduces the asymptotics of det M N (&) to the asymptotics of 
the (skewsymmetric) Toeplitz determinant det T 2 jv(aer). In order to analyze this 
Toeplitz determinant we cannot rely on the (original) Fisher-Hartwig conjecture 
because it breaks down in this case. However, in Section || we will proof a limit 
theorem saying that the quotient 



converges - under certain conditions on a - to a nonzero constant. In order 
to prove this limit theorem we make heavy use of the machinery that has been 
developed in || in order to prove the Fisher-Hartwig conjecture. 

Thus, up to this point, we have reduced the asymptotics of det Mjv(a) to the 
asymptotics of det T 2 Ar( a )- The Toeplitz determinant detT 2 Ar(a) is (generically) 
of a kind for which the Fisher-Hartwig conjecture holds. In Section [| we will 
therefore recall the Fisher-Hartwig conjecture in the form as it has been proved 
in |9|. Moreover, we specialize it to the distributions (namely, even distributions 
of Fisher-Hartwig type) that we are interested in. In Section || we combine all 
the previous results and obtain the asymptotics of det Mat (a) for even distri- 
butions a of Fisher-Hartwig type, which satisfy appropriate conditions on the 
parameters. 

In Section ^ we finally specialize the quite general result of Section [6] to even 
piecewise continuous functions. We thus obtain the asymptotics of det M^^) 
for a certain class of even piecewise continuous functions. This result together 
with the results that are known from [Q suggest a conjecture about the asymp- 
totics of detMjv(0) for quite general piecewise continuous, not just those nec- 
essarily even. We end the paper with the conjecture. 



(det Mat (a)) 2 =detT 2JV acr). 



(4) 



det T 2 jv(acr) 
detT 2A r(a) 



(5) 
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2 Known results for determinants of Toeplitz 
plus Hankel matrices 

Let us begin by recalling some of the results concerning the asymptotics of the 
determinants of the matrices M/v(</>) that have already been established in Q|. 

We first consider the case of continuous and sufficiently smooth generating 
functions </>, where an analogue to the Strong Szego-Widom Limit Theorem 
holds. In order to be more specific about the smoothness condition, let us 
consider the Besov class B\, which is by definition the set of all functions b £ 
L X (T) such that 

||&||m := / \\ \b{e lx+ly ) + b{e lx - ly ) -2b{e lx )\ dxdy < oo. (6) 



— 7T V J — TT 

It is known that B\ forms a Banach algebra with the above norm and is con- 
tinuously embedded into the Banach algebra of all continuous functions on T. 
By G\B\ we denote the set of all nonvanishing functions in B\ with winding 
number zero. The set G\B\ can also be characterized as the set of all functions 
b which possess a logarithm log b in B\ . 
For b £ G\B\, the constants 

G[&] = exp([log6] ), (7) 

E[b] - exp f^frpog&lfcpog&U j , (8) 



(oo 1 OO \ 

$>g6] 2fe - 1 --5>[log&] 2 fc , 
fe=l k=l I 



(9) 



are well defined, where [log6] n stand for the Fourier coefficients of log b £ B\. 
Moreover, for b £ G\B\, the functions b + , b- £ GiBl are well-defined by 

b±(t) = exp (jtt ±n l\ogb] ± }j , tel (10) 

Note that b(t) = b-(t)G[b]b+(t), t £ T, is the normalized canonical Wiener-Hopf 
factorization of the function b. 

The analogue to the Strong Szego-Widom Limit Theorem for the determi- 
nants det Mat (</>), which has been established in Corollary 2.6], now says that 
if b £ G X B\, then 

det M N (b) - G[b] N E[b]F[b] as N -> oo. (11) 
In the case of even functions b £ G\B\ this simplifies to 

det M N (6) - G[b] N E[b] as N -> oo, (12) 
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where E[b] is the constant 

(-.00 00 \ 

-Y / ^ogb} 2 k + ^°S b hk-i) ■ (13) 
fe=l k=l ) 

In order to discuss the asymptotics for the case of piecewise continuous 
generating functions <j), let us introduce the functions 

t/3,e (e i9 ) = e^ e - e °-*\ < 9 - 8 < 2tt, (14) 

where (3 G C and 8q G (—it, it]. The piecewise continuous functions that we 
consider arc of the form 



(e ie ) 



b(e i6 )l[t^e r (e i6 ), (15) 



r=l 



where Q\,...,Qr G (— 7T, 7r] are distinct numbers determining the location of the 
jump discontinuities and are complex parameters determining the 

"size" of the jumps. The function b is usually assumed to belong to G\B\. 

As is known from the theory of Toeplitz determinants, a key ingredient 
for the determination of the asymptotics of the determinants are localization 
theorems (see, e.g., ||). A localization theorem for the determinants det Mjv(^) 
with </> being piecewise continuous has been established in @, Theorem 5.11]. 
This localization theorem reduces the asymptotics of det Mm{4>) for "general" 
piecewise continuous functions ( |l5| ) to the asymptotics for particular piecewise 
continuous functions. 

Theorem 2.1 (Localization Theorem) Let 4> be a function of the form 

R 

0(e* e ) = b(e w )Me l9 )4-(e l9 ) Y[ <Me l9 ), (16) 

r=l 

where b G G\B\, 4> + = i/3 + ,o, 0- = t(3-,n an d 4>r — tp+ g tg- f or I < r < R. 
Suppose that 6\, . . . , Or G (0, it) are distinct numbers and that /3±, fif, . . . , ftf; G 
C are such that 

(a) -1/2 < Re/3+ < 1/4 and -1/4 < Re /3_ < 1/2, 

(b) |Re/3+ < 1/2 and |Re#r| < 1/2 and |Re(/3+ +/?")! < 1/2 /or eac/i 
1< r < R. 



The 



,. det Mat (0) 

hm = ff , 

jV--ao -R 

det M N (b) det M Ar (^+) det M N (4>-) J| det M N (((> r ) 

r=l 
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where 

H = & + (l)^+6_(l)-^+6 + (-l) 2 / 3 -6_(-l)-^-2 3 ^^ 

n6 + (t r )^+^6_(t f .)-^6 + (tr 1 ) /s ' ++ ^6-(*r 1 )- /,; 



H 

X 
r=l 

-R 

X 

r=l 

.R 

X 

r=l 



JJ(1 - t r )/M^ + + 2 ^)(l - t-l)/3 + (2/3 r + +/3D 
r=l 

^(1 + tr f-{0t+Wr){l + t -y-(20++Pr) 
--1 

H (i - trt a )KK + rtK+Krt(i - t-H^frPi+^K+^Pi 

r<S<R 



X 

l<r<s<fl 



X 

l<r<s<fi 



ifere i r = e l6,r , 1 < r < -R, and &± are i/ie functions (|10|). 

The asymptotic behavior of detM/v(</>) with the generating function + = 
i/3 + ,o an d = tp_ jV) respectively, has also been determined (see Theorem 6.2 
and Theorem 6.3 in [0]). These functions have one single jump discontinuity at 
the points 1 and — 1, respectively. 

Theorem 2.2 Let j3 G C \ Z. Then 

.... detM N (t , o ) a 3^G(|-/3)G(l-/3)G(l + /3) 

(a) hm , = (27r)2 2 2 — = , 

1 ' N -^-§ { ' G(±) 

(b) Um detM^) _ {ti ^ 2 s£G(l-0)G(l-0)G(l + P) 



2 ' 



In these asymptotic formulas the Barnes G-function G(z) appears [g, 
which is an entire function defined by 

00 / 1 \ 

G(l + z) = (27r)te-™- CE ^ J] Ml + I) e -*+£j (17) 



fc=i 

with G^ being Euler's constant. 

In the case of the generating functions <p r — t^+ e t^- _ g , which have two 
jump discontinuities at a point of the unit circle and its complex conjugate, 
the asymptotic behavior is only known in particular cases. One case is that 
where either /?+ = or f3~ = 0, i.e., where the function has exactly one jump 
discontinuity at a point in T\{1, — 1}. Here the result is taken from |4|, Theorem 
4.5]. 
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Theorem 2.3 Let 9 G (-tt, 0) U (0, tt) and f3 G C be such that \Ref3\ < 1/2. 
Put t = e l(>0 . Then 



det M N (tg fio ) 



2.1 



to the 



jt'oo ^"^;r o ' = g(i - p)G(\ +0) (i - 1^ 1 ) 2 + ^ (i + v 1 ) 2 * . 

Finally, another, even more particular case of a function with two jump 
discontinuities at i and —i and the same size of the jumps has been established 
[|] if one combines Theorem 7.4 and Theorem 7.5. 

Theorem 2.4 Let j3 € C \ Z. Tften 

det M wit a 2Ltg _n) .„ 2 „ 
Jim , 3 ffl ^ = 2 4/3 G l - 2/3 G 1 + pf . 

N^oo iV 

For later use, let us specialize the localization theorem (Theorem 
case of even functions <f> which are of the form (|l6|). 

Corollary 2.5 Let 4> be a junction of the form 

R 

<t>(e ie ) = b(e ie )l[Me i9 ), (18) 

r=l 

where b S Gi_B| is an even function and (f> r — ta r e r ^—dr —Or f or 1 < r < i?. 
Suppose that 9i, . . . ,9r G (0, 7t) are distinct numbers and that 0i, . . . , /3r G C 
are smc/i i/iai |Re/3 r | < 1/2 /or eac/i 1 < r < i?. Then 

hm detMjv( ^ = H, 

N- 



det M N (b)Y[ det M N (cf) r ) 



r=l 



where 



-Pr 

/ -\- \ u r ) L/ — J 
r=l 

x n (i-trt.r^a-^V)-*-* 

1<?-<s<,R 

x n (i-trtr^-a-^ 1 *.)^. 

l<r<s<7?. 

iJere t r = e* er , 1 < r < R, and b± are the functions (|1 



Proof. We apply Theorem 2.1 with the parameters j3± = and (3^ = ±/3 r . We 
remark also that o + = 6_. □ 
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3 Preliminary results for determinants of 
symmetric Toeplitz plus Hankel matrices 

The first step in order to determine the asymptotics of the determinants of sym- 
metric Toeplitz + Hankel matrices M/v(</>) is to express these determinants by 
means of determinants of skewsymmetric Toeplitz matrices. The identity as it 
appears in the following theorem has been stated explicitly in JlJ, Lemma 18], 
but is already implicitly contained in |n], Lemma 1] and jl4| Proof of Theo- 
rem 7.1(a)], where it has been proved. A different, self-contained proof has been 
given by the authors in ||. 

Theorem 3.1 Let {a^^L-oo be a sequence of complex numbers such that 
a_„ = a n . Let 

n 

c n = ^2 a k for n>0, (19) 

fc=-n+l 

and put cq — and c_„ = —c n . Then 

(det (aj-k + a j+k+ i )f fc= 1 Q ) = det (cj- k )flzl ■ (20) 

The matrices appearing in ( |2p| ) are a symmetric Toeplitz + Hankel matrix 
of the same kind as (|l|) and a skewsymmetric Toeplitz matrix. If we are trying 
to rewrite this identity by using the standard notation (||) for Toeplitz and 
Hankel matrices where the sequences {a n } and {c n } are the Fourier coefficients 
of functions n,c£ L 1 (T), we face the difficulty that this is in general not possible. 
Consider for instance the simplest case where a(t) = 1, i.e., ao = 1 and a n — 
if n 7^ 0. Then we obtain c„ = sign(n), and obviously, there does not exist a 
function c 6 L 1 (T) with such Fourier coefficients. 

A way out of this situation is to consider distributions on the unit circle in 
place of functions in L : (T) and to take their Fourier coefficients as the entries of 
the Toeplitz and Hankel matrices. For this purpose we need several preliminary 
results. Apart from basic issues, the following has all been stated in jl0| and 
proved in M. 

Let V = C°° (T) be the linear topological space of all infinitely differentiable 
functions defined on the unit circle. By V we denote the set of all distributions 
on the unit circle, i.e., linear and continuous functionals on T>. The Fourier 
coefficients of a distribution a£D are defined as 

a„ = a(x-„), (21) 

where \n S L) is the function Xn(t) = t n . There is a natural identification 
of functions a £ L X (T) with a subset of distributions. It is established by the 
mapping a 6 L 1 (T) naeP' where 

a(/) = ^ 1^ a ^ l8 )f^ e ) d0 > f eV - ( 22 ) 
2tt Jo 
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This definition ensures that the Fourier coefficients of a and a are the same. We 
also remark that there is a one-to-one correspondence between V and the set 
of all at most polynomially increasing sequences {fln}^.^, which is given by 
associating to a G V the series {a n }^' = _ 00 of its Fourier coefficients. 

Let a G T> and h E V . Then the product of a and b is the distribution 
ah G V which is defined by 

(ab)(f)=b(af), feV. (23) 

Let K be a compact subset of T. We denote by C°°(T \ AT) the set of 
all infinitely differentiable functions on T \ K. By CJf(T) we refer to the set 
of all functions / G T> which vanish on an open neighborhood of K . The 
product of a function f G C°° (T \ AT) with a function g G C^? (T) is a function 
fg G C% (T) C V by putting (fg) (t) = for t G K . 

We will proceed with some definitions that are not so quite common, but 
necessary for our considerations. They are taken from g, |lQ| . Let V'(K) stand 
for the set of all distributions a for which there exists a function a G C°° (T \ K) 
such that 

a(/) = i a(e ie )f(e ie )d6 (24) 
27r Jo 

for all / S C^(T). The function a s C°°(T\K) is uniquely determined by a and 



called the smooth part of the distribution a. Definition (24) can be rephrased 
by saying that /a = fa for all / S Cjj? (T), where the left hand side of this 
equation is a distribution £>' and the right hand side is a function in T>, which 
are identified in the sense of (]22j). 

Next we are going to show that one can define - under certain assumptions 
- the product of two distributions. Let M and N be compact and disjoint 
subsets of the unit circle. Given a G V(M) and b G T>'(N) with smooth parts 
a G C°°(T \ A/) and 6 G C°°(T \ AT), let ab G P'(M U AT) be defined as 

ab = (6/ b )a + (a/ Q )b, (25) 

where f a G C^(T) and f b G C^(T) are such that / a + / 6 = 1. This definition is 
independent of the particular choice of f a and fb- Moreover, ab has the smooth 
part abe C°°(T\(M UN)). 

Given a function a defined on (a subset of) the unit circle, we define the 
function a by a(t) = a(t _1 ), t e T. In accordance with this definition, given a 
distribution a G V , let a G V stand for the distribution with Fourier coefficients 
a_„ = a„. A distribution a will be called even (odd) if a = ±a. A function 
a will be called even (odd) if a = ±a. Finally, if AT is a subset of T, put 
K = {t G T : t^ 1 G AT}. If K = AT, we call K a symmetric subset of the unit 
circle. 



In the reformulation of Theorem 3.1, the following distribution will play a 



role. Let a G V be the distribution which has the Fourier coefficients 

cr n = sign(n). (26) 



9 



Moreover, let a G C°°(T \ {1}) be the function 

"(*) = ^lf (27) 
Remark that both the distribution a and the function a are odd. 
Proposition 3.2 The distribution a is in 2?'({1}) and has the smooth part a. 
Proof. For / G C { °J } (T) we can write f(t) = (1 - t)h(t) where heV. Then 

[/c]n = ^ fn-k — ^ /n-fc 
fe>0 fe<0 

= / — h>n-k-l) — ~^^{h n -k ~ h n -k-l) 

k>0 fc<0 

= /l„_i + /l„ = [/l(t)(l + i)]n = [Mi- 

Note that /ijt converges to zero sufficiently fast. This completes the proof. □ 

Given a distribution aeD' with Fourier coefficients {a„}^L_ OOJ define the 
N x N Toeplitz and Hankel matrices by 

7V(a) = (a,-*)^,, H N ( a ) = ( aj+k+1 )^ Q . (28) 

This definition is in accordance with (|J). Moreover define 

M N (a) = T N (a) + H N (a). (29) 



Now we are ready to give the desired reformulation of Theorem 3.1 



Theorem 3.3 Let K be a compact and symmetric subset o/T\{l} 7 and assume 
that a G T>'(K) is an even distribution. Then 

detT 2A r(aer) = (dct Afjv(a)) 2 . (30) 

In the above c = acr G T>' {K U {1}) is an odd distribution. 

Proof. Since K and {1} are disjoint compact sets, the distribution c = acr is 
well defined. In the definition c = (cr/ CT )a+ (af a )o~ we may assume without loss 
of generality that f a and f a are even functions. From this it follows easily that 
cr/cr is odd and af a is even; thus both (af a )a and (a/ a )cr are odd. Hence c is 
odd. 

Next write f a {t) = ^(l-tXl-t- 1 ). Then (of a ){t) = {l+t){l-t- x )g{t) = 
(t - We obtain that 

[(o-/<r)a] n = [((* - t- 1 )g)a] n = [ga] n _i - [c/a]„ + i. 
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Moreover, using the fact that f a {t) = g{t)(l — t)(l — t x ) and keeping track of 
the cancellation, it follows 

n n 

if* a h = (-[9^-1 + 2[ga] fc - [ga] k+ i 

k— — n-\-l A;— — n+1 

= -[ga} n + [ga]„_i + [ga]- n - [#a]_ n _i 
= [ga]n-i - [ga] n +\- 

Here we have also used that ga is even. From these two identities we obtain 

n 

[(aU) a ] n = t-M*- ( 31 ) 

fc=— n+1 

On the other hand, since af a = a/ a is even, 

n 

[(afa)cr)] n = y^[afg] n -k - ^[a/o]n-fc = Yj [ a /J fc - ( 32 ) 
fc>0 fc<0 fc=— n+1 

Combining @ and © yields 

n 

fe=-n+l 



Together with Theorem 3.1 this completes the proof. □ 
Finally we will need the following result. 

Proposition 3.4 For each N > 1 we have detT 2 jv(c) = 1 and 

T^{a)=T 2N {v) (33) 

where v is the distribution with Fourier coefficients v n = sign(n)(— 1)™. 



Proof. By Theorem 3.3 we have detT 2 N(a) = (detMjv(l)) , where obviously 



Mjv(1) = In- The formula for the inverse of T 2N (a) can be easily checked. □ 

4 A limit theorem for determinants of skewsym- 
metric Toeplitz matrices 



Theorem 3.3 reduces the computation of the asymptotics of detMjv(a), which 
is what we are interested in, for certain distributions a, to the computation of 
the asymptotics of det T 2 m (ca) ■ 

At first glance one might think that the asymptotics of det T2 at (era) could 
be obtained from the predictions of the Fisher-Hartwig conjecture, which was 
proved in H. Unfortunately, since T 2 N(aa) is a skewsymmetric Toeplitz matrix, 



II 



the Toeplitz determinant belongs to those classes of functions where the Fishcr- 
Hartwig conjecture breaks down. It might be the case that the asymptotic 
behavior fits with the still unproved generalized conjecture |?], [T(|. However, 
distributions of the kind era appear here probably for the first time in connection 
with Toeplitz determinants and since previous results do not include this setting, 
previous techniques must be modified. 

To that end, the purpose of this section is to prove that under certain as- 
sumptions on the distribution a the expression 

detr 2Af (acr) 

detT 2Ar (a) 1 ' 

converges to a certain (explicitly given) nonzero constant as N — > oo. Although 
we cannot rely on the main results of Q (see also flOfl), i.e., the Fisher-Hartwig 
conjecture, we will very heavily rely on the machinery and several auxiliary 
results established in || . 

Let us proceed with recalling the necessary definitions. For fj, G R let £ 2 
stand for the Hilbert space of all sequences {x n ]'^ =0 of complex numbers for 
which 




KK?=oL:= Wl + nHa^ < oo. (35) 



For fii > // 2 the space is continuously and densely embedded in i 2 . 

The Toeplitz and the Hankel operator generated by a G V are the one-sided 
infinite matrices 

T(a) = (a,- fe )^=o , H(a) = (a,- +fe+1 )~ =0 , (36) 

where a n are the Fourier coefficients of the distribution a. For each acD' there 
exist a (sufficiently large) fii and a (sufficiently small) /i 2 such that the matrices 
T(a) and H(a) represent linear bounded operators acting from into £^ 2 . 

The situation of the boundedness of Toeplitz and Hankel operators gener- 
ated by functions in T> was established in the following lemma taken from 
Sect. 6.2]. 

Lemma 4.1 For each fi, \x\,\xi S K and a 6 T>, the operator T(a) is bounded 
on £(£^,^) and the operator H(a) is bounded on £(£^,1^). 

We define the following finite rank operators acting on I 2 -. 

P N : (x 0l xi,x 2 , ...)>-* (x ,Xi, . . .,x N -2,Xn-i,0,0, . . .), (37) 
W N : (x Ql xi,x 2 , ...)>-> (xfj-x,x N -2, . . . ,xi,a;o,0,0, . . .). (38) 

Obviously, — — Pn and WnPn = PnWn = Wn- If we consider the 
matrix 2V (a) as acting on the image of the projection Pjy i n the space then 
T N (a) = P N T{a)P N . Moreover, 

W N T N (a)W N = T N (jk). (39) 
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Recall that a is the distribution with the Fourier coefficients a„ = a_ n . 

For our purposes we need to single out two additional classes of distributions. 
Let T>' + (T>'_, resp.) stand for the set of all distributions a£D' for which a„ = 
for all n < (n > 0, resp.). These two sets form commutative algebras with 
a unit element e(t) = 1. For a, b <E D' + , the product c = ab is defined by 
stipulating c n — for n < and 

n 

c„ = ^2 a »-fe b fc for n>0. (40) 

fe=0 

For a, b G T>'_ , the product c = ab is defined by stipulating c n = for n > 
and 

o 

c„ = ^ a„-feb fc for n < 0. (41) 

k—n 

This definition of a multiplication is compatible with that of ( p5[ ) whenever both 
are defined. Let QT>'± stand for the group of all invertible distributions in T>±. 
Moreover, we put V'±(K) = V'±n V'{K) and let gV ± (K) stand for the group 
of all invertible elements in T>'±(K). One can show that if a G QV' ± {K) has the 
smooth part a, then a is an invertible element of C°°(T \ K). 

There are some obvious relations between the distributions a and a. For 
instance, if a G QT>' + (K), then a G QV'_(K). Moreover, if a has the smooth 
part a, then a has the smooth part a. 

Let Hi and H 2 be Hilbert spaces. We consider sequences {CW}^ =1 the 
elements of which are well defined linear bounded operators (or matrices) Cn '■ 
Hi — > Hi for all sufficiently large N. Let O(g) with g G R stand for the set of 
all such sequences for which 

\\CN\\c{H 1 ,H 2) =0{N e ) asiV^rc. (42) 

The dependence of O(g) on Hi and H^ will not be displayed in the notation. 
We also use the notation O(g) in order to denote any sequence of this type. In 
this sense, C N = C + O(g) means that {C N - C}^ =1 G O(g). 

Now let Hi, H2, Hi and Hi be Hilbert spaces and go, gi, gi G R. We denote 
by O(go, Qx, Q2) the set of all sequences {Cn}^=i of 2 x 2 block operators for 
which 

C "=(?n( 0O \?n( gi \) ■ Hi®Hi^H 2 ®%. (43) 

We will also use the notations 02(g) and Oi(g). The only difference in 
comparison with O(g) is that we consider the convergence in (^|) in the Hilbert- 
Schmidt and in the trace class norm, respectively. Likewise, we will use the 
notations 02 (go, Qi, Q2) and Oi(g , QuQ2)- 
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Given a € V , assume that Tjv(a) is invertible for all sufficiently large N, 
and introduce the following sequences of operators of 2 x 2 block form: 



HRH N (a) 



H(a)P N T^(a)P N H(a) H(a)P N T^(a)W N H(a) 

H(a)P N T N 1 (a)W N H(a) H(a)P N T^(a)P N H(a) 

T(a) H(a X - N ) 
H(ax-N) T(a) 



(47) 



Here, as before, x~N(t) = t~ N , t G T. These sequences of operators are consid- 
ered from ^ © £^ into © with /xi sufficiently large and /12 sufficiently 
small, which ensures the boundedness of the operators. 

Now we are prepared to define the notion of ^-convergence, which has been 
introduced in || Let a 6 V , a + G QT)' + and a G QT>'_ be distributions, 
and let go, g\, g 2 G K. Wc say that the distribution a effects 1Z- convergence with 
respect to [a + ,a_] and (go, Qi, g 2 ) if there exist fj,% > and \x 2 < such that 



fljv(a) = diag [T(a; 1 )T(a: 1 ), ^(al^T^; 1 ) j + O(g , Ql , g 2 ), (48) 

RH N (a) = diag ^(a; 1 )// (a+), T^a: 1 )^^)) + O(g , Qi, Qi), (49) 

HR N (a) = diag (ff (a_)T(aI 1 ), H(a+)T(a; 1 )) + O(g , Qi, Q 2 ), (50) 

HRH N (a) = - diag (V(a_)T(a+), T(a+)T(a_)) + O(g , gi, g 2 ), (51) 



where these sequences are considered from © £f, into £?. © . 

fji\ f-L\ fj,2 l- L 2 

We will need the concept of "^.-convergence for a distribution a which is even. 
This particular case gives rise to some simplifications. For any distribution 
a G T>' , the following statement are equivalent: 

(1) a effects ^-convergence with respect to [a+,a_] and (go, g\, £2); 

(2) a effects 7?.- convergence with respect to [a_, a+] and (go, g 2 , gi)- 

In fact, in order to prove this equivalence, one has only to pass to the transpose 
in equations (f48|)-(|5l|). Hence for an even distribution a we can replace g\ and 
02 by minjpi, 02}, i-e., we may assume that g\ = g 2 =: g. Moreover, since the 
distributions a + and a_ are uniquely determined up to a nonzero multiplicative 
constant, we may assume without loss of generality that a = a + . 
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This last remark gives some motivation for the assumptions in the following 
theorem. In this theorem we establish the asymptotic formula for (|54j). 

Theorem 4.2 (Limit Theorem) Let K be a symmetric and compact subset 
o/T\{l, — 1}. Moreover, assume that 

(i) a £ V (K) is an even distribution with the smooth part a £ C°°(T \ K); 

(ii) a + G QT>' + (K) is a distribution with the smooth part a + £ C°°(T \ K); 

(iii) a{t) = a + {t)a + {t) for all t £ T \ K; 

(iv) go < and g < 0; 

(v) a effects 1Z- convergence with respect to [a+,a+] and (go, g, g). 
Then 

detT 2N (aa) o+(l) 



detT 2N (a) o+(-l) 



0(iV max ^o,e}) as 7v -> oo. (52) 



The rest of this section is devoted to the proof of this theorem. Once again 
we need to quote some auxiliary results. 
Let us introduce the notation 

*£=&©C * W =C */*=C (53) 

This notation is convenient in the sense that it reflects the condition that [i\ 
is a sufficiently large and \x 2 a sufficiently small real number, which we will 
encounter in what follows. 

The following proposition is taken from || Proposition 8.1] with a slight 
change of notation. As before, we will denote a distribution by a bold letter and 
its smooth part by the same non-bold letter without mentioning it explicitly. 

Proposition 4.3 Let K be a compact subset of T, and let a £ D'(K) and 
f £ CjJ? (T). Then for all sufficiently large fii > and all sufficiently small 
ji 2 < the linear operators 

#i(a, /) = (ff (a), (o/) - T(a)ff(/) ) , (54) 

W,a)=(^ a )- ( ^>W) (55) 

are bounded on the spaces 

H X (a, /) : X£ - X^ , ff 2 (/, a) : -> X£ . (56) 

In order to provide some meaning to these operators we remark that with ap- 
propriately chosen [i\ and fi 2 , these operators can be embedded into the spaces 

Hi (a, /) : X^ -> X M2 , H 2 (/, a) : X^ -> X£| . (57) 
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In this case, these operators can be written as a product: 

H 1 ( a ,f)=H( a )(l,T(f)), (58) 



ff 2 (/,a) = ( T{ P )H(S). (r,!)i 



I 

The next result is taken from [|9|, Corollary 8.8]. 

Proposition 4.4 Le£ K\ and K 2 be disjoint and compact subsets off, and let 
&i e V'(Ki) and /* G Cf.(T) (i = 1,2) 6e suc/i i/iai /i + / 2 = 1. 77ien 

Tjv(aia 2 ) = T^(a 1 )T iV (a 2 ) + Pjvffifai, fi)H 2 (f 2 , a 2 )Pv 

+ WVffi(ai,/ 1 )iT 2 (/ 2 ,a 2 )W i v. (60) 

We remark that the linear operators occurring in ( |60| ) are bounded on appro- 
priately chosen spaces. Moreover, ( |60| ) represents a generalization of the well- 
known identity 

T N {a ia2 ) = T N ( ai )T N {a 2 ) + P N H{ ai )H{h 2 )P N + W N H{a x )H{a 2 )W N 

due to Widom [ jig , which holds for functions oi, a 2 € L°°(T). 

Next, we consider the functions £i(i) = 1 — t^ 1 and £-i(t) = 1 + t — 1 . 
These functions can be identified with distributions £i G T>'_ and £_i G 
respectively, in the sense of (p2|). Obviously, both £i and £_i belong to QT>'_. In 
fact, the inverse distributions ^j -1 and £Zi are given as follows by their Fourier 
coefficients: 

^ z i]n = {(-i)« ifn<o. (62) 



Proposition 4.5 TTie following statements are true: 

(a) £i G ^I>L({l})i ond C? 1 ^ e smooth part ^j^ 1 ; 

f&) £_i G QTy_(\— 1}), and /ias i/ie smooth part 

Proof. Since £i = £i in the sense of (p2[), we even have £i G 2?'(0) and £i 
has the smooth part £i. It remains to show that ^f 1 is contained in V'({1}) 
and has the smooth part £j~ . Indeed, let / G C^j(T). Then we can write 

/ = £i5 where 5 G C°°(T). It can be checked easily that £f £1 = 1. Hence 
£f / = g = Cf 1 /- This completes the proof of part (a). Part (b) can be proved 
analogously. □ 

The following corollary, in which we define another distribution h, is a simple 
consequence of the previous proposition. 
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Corollary 4.6 The distribution h = ^-j - 1 ^ i is contained in QT>'_{{— 1, 1}) and 

has the smooth part h — ££~ (,—x- The inverse distribution hr 1 equals £i£Zi and 
has the smooth part h^ 1 — £i£Zi- 

It is easy to compute the Fourier coefficients of h and h _1 : 

if n > 

[h]„ = { 1 if n = (63) 
2 if n < 0, 

if n > 

1 if rc = (64) 
2(-l)™ if n < 0. 

There is (for our purposes) an important relation between operators con- 
taining the distributions a and h, which is given in the following proposition. 

Proposition 4.7 Let / G C^°_ 1 1} (T). TTiera 

H 2 (f,a)P 2N T-i(a) = H 2 (f,a)W 2N T^(a) = 

1„ „ < 1 

2 

Proof. We first prove that 



= -H 2 (f,h)P 2N T-i(h) = --H 2 (f,h)P 2N T-i(h)W 2N . (65) 



H(a)P 2N T^(a) = H{a)W 2N T^(a) = 

= ~H(h)P 2N T-J(h) = -~H(h)P 2N T-i(h)W 2N . (66) 

In this identity we do not need to worry about the boundedness on certain 
spaces since both the left and right hand side is an infinite Hankcl matrix times 
a finite rank matrix. Let x = (1,1,.. .) T denote an infinite column vector and 
x 2 N — (1, ■ ■ ■ , 1) T a finite column vector of size 2N. Then 

H(a)P 2N = H{d)W 2N = -xxl N , H(h)P 2N = 2xx T 2N . (67) 



A moments thought shows that (|66J) is proved as soon as 

- x 2n^ 2 n W = xInT*& (h) = ~x^ N T^ (h)W 2N (68) 

is established. However, this is just a straightforward calculation. We have to 
observe that T 2 ~^(h) = T 2 n()i~ 1 ) with the Fouri er co efficients given by ( |64|) and 
moreover that 

(a) — T 2 n{v) by Proposition p. 4j 
Having proved (|6^), we take into account the identity (59) and formula ( |65|) 
follows by a density argument of the Hilbert spaces under consideration. □ 

In §, Formula (8.64)] , the following 2x2 block operators acting on Xfe 
with sufficiently large \i\ and small \i 2 were defined: 

HSH N (f 2 , a2 , ai> h) (69) 
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Moreover, it has been shown |)| Formula (8.139)] that 

HSH N (f 2 , a2 , ai ,f 1 ) = HZ N (f 2 , a2 J 2 )YH N (f 1 , ai ,f 1 ), (70) 

where HZn(. . .) and YHn(. . .) are linear bounded 2x2 block operators acting 
on © X^ , which we are not going to define here. The important result 
concerning these operators is the following asymptotic formula, which is taken 
from H Proposition 10.4 and Proposition 10.5]. 

Proposition 4.8 Let K be a compact subset of T, / G C^?(T), and assume 
that a G T>'(K), a± G QT>'±{K) such that a = a + a_ holds for their smooth 
parts. If a effects 1Z- convergence with respect to [a+,a_] and (go, Qi, Q 2 ) G K 3 , 
then 



y j ff JV (/,a,/) = diag \TH(f,a+J), TH(f, a ^J))+O 2 (g ,gx,Q 2 ), 
HZ N (f,aJ) =diag fflT(/,a_,/), HT(f, a+, /)) + O 2 (q , Qx, Q 2 ). 



In the previous proposition linear operators TH{. . .) and HT(. . .) appear, 
which were defined in || Formulas (9.11) and (9.12)]. These operators are 
Hilbert-Schmidt operators on the space X^ for all sufficiently large [i\ > and 
all sufficiently small \i 2 < (see ||, Proposition 9.2]). 

Moreover, the asymptotic operator relation stated in the previous propo- 
sition has to be understood in the way that the operators act on the space 
X^ © Xfi* , where fix > is fixed and sufficiently large and fi 2 < is fixed and 
sufficiently small. 

Now we are prepared to give the proof of Theorem 4.2. 



Proof of Theorem 4.2. We start from Proposition 4.4 with ai = a, K\ = K, 
a.2 = cr, K 2 = { — 1,1}. Since K and {1,-1} are symmetric subsets of K, we 
can assume without loss of generality that f\ and f 2 are even functions. Then 

T 2N ( a <j) = T 2N { a )T 2N {a) + P 2N H x { a , ft)H 2 (f 2 , a)P 2N 

-W 2N Hx{ a ,fx)H 2 (f 2 ,a)W 2N , (71) 



where we have also used that a is even and a is odd. By Proposition 3.4 the 
inverses of T 2 n{o~) exist for all N. Since the distribution a effects 1Z- convergence, 
the inverses of T 2 n(sl) exist for all sufficiently large N. Hence 

T-i( a )T 2N ( a a)T^(a) = P 2N + T 2W 1 (a)P 2A ri/ 1 (a, h)H 2 {f 2 , a)P 2N T^(a) 

- (a) W 2Jv ffi (a, fx)H 2 (f 2 , a)W 2N T^ (a) . 



From Proposition 4.7 it now follows that 
T^( a )T 2N ( a a)T 2 -i(a) = P 2N + ~T^( a )P 2N Hx( a , fx)H 2 (f 2 , h)P 2N T^(h) 

- l -T 2 ^ ( a )W 2N Hx (a, fx)H 2 (f 2 , h)P 2N T 21 $ (h) . 



18 



Taking determinants, observing that det T 2 n(<j) = 1 by Proposition 3.4 and 
using the formula det(I + AB) = det(I + BA) for determinants, we obtain 

detT 2A r(a ( j) 



detr 2JV (a) 



1, 



det \P 2N + -T-A( a )(P 2N - Wj^ff^a, fi)H 2 (f 2 , h)P 2N T~A(h) 
det fl+^H 2 (f 2 ,h)P 2N T 21 ^(h)T-^( a )(P 2N - W 2N )H 1 {ai,f 1 : 



det ( /+ -H 2 {f 2 ,h)P 2N T-i(h){P 2N - W 2N )T- 1 (a)P 2N H 1 {aJ 1 [ 



1 

r 

Here we have also used formula (|39|) and the fact that a is even. Again from 
Proposition 



4.7 it follows that 



detT 2j y(aa) 
detT 2JV (a) 

= det (/ + H 2 (f 2 , h)P 2 ^T 2 ] v 1 (h)T 2 ^(a)P2 J v-ffi(a, fx)) 

= det 0)HSH 2N (f 2 ,h,aJ 1 ) (Ijy (72) 

Note that in the last formula the (1, l)-block entry of HSH 2 iy(f 2 , h, a, /i) ap- 
pears (see formula (^9|)). 

From the assumption on a and - as concerns the distribution h - from 
Proposition |4.5|, C orollary 4.6 and [|[ Theorem 13.1] (see also the remark made 
after Theorem |5.2| below) we know that 

(1) a effects ^-convergence with respect to [a+, a + ] and (go, g, g); 

(2) h effects ^-convergence with respect to [1, h] and (ui, 0, ui) for each wel. 
Hence, by Proposition [D^ we can conclude that 

'TH(f 1 , a+ J 1 ) + O 2 (g ) 2 (g) 

(HT(f 2 ,hJ 2 ) + 2 (Lu) O a (0) 
^ 2 (u) HT(f 2 ,lJ 2 )+0 2 (u>) 

This in connection with ( |70| ) yields 

I 




YH 2N {h,a,h) 
HZ 2N (h,hJ 2 ) 



(/, 0^HSH 2N (f 2 ,h,a,f 1 ) 

= (HT(f 2 ,h,f 2 ) + O 2 (cj))(TH(f 1 , a+ J 1 ) + O 2 (g )) +O 2 (0)O 2 (g). 



Noting that the operators HT(. . .) and TH(. . .) are Hilbert-Schmidt and choos- 
ing lo sufficiently small, this implies 

(/, Q)HSH 2N {f 2 ,\x,aJ x ) (I) 

= HT(f 2 , h, f 2 )TH(f 1 ,a+,f 1 ) + Oi(max{eo, *?}). (73) 
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Formulas ([72|) and (|73|) give 

<l0t T ' 2X Uln ~ ! dct (I + HT(f 2 , h, f 2 )TH(f ua+ , A)) + 0(i\r»<*"«>). 



detT 2Ar (a) 

It remains to show that the above operator determinant equals the constant 
a+(l)/a+(— 1). From || Proposition 9.10(b)] we obtain 

det(/+ffT(/ 2 ,h,/ 2 )Tff(/ 1 ,a + ,/ 1 )) = lim E(h r ei + , h r h) (74) 

r— — 

where /i r a + and /i r h (0 < r < 1) are the harmonic extensions of the distributions 
a + and h, i.e., 

OO CO 

(h r a + )(t) = r n t n [ a+ ] n , (h r h)(t) = r"r"h_„, (75) 

n=0 n=0 

and E{. . .) is the constant defined by 

E(<j> + ,4>_) = cxp ^n[log^ + ]„[log0_]_ n J . (76) 

In this connection we remark that the harmonic extensions of distributions 
contained in QV' + or in QT>'_ are always functions in T> = C°° (T) , which possess 
a continuous logarithm on T for each < r < 1. Indeed, if b € QT> ±, then the 
harmonic extensions are multiplicative, and consequently (h r h)(h r b ) = 1. 
The harmonic extensions depend uniformly on r and are constants for r = 0. 
Hence the functions h r h are nonzero on all of T and have winding number zero. 
In order to compute E(h r a. + 1 h r h), observe first that 



Hence 



[log(/ lr h)]_„ = - - ^-i-, n > 1. (77) 
n n 

We obtain that E(h r a. + , h r h) is the exponential of 

OO 

J2 ([logh r a + } n r n - [log/i r a + ]„(-r)")) 

n=l 

= (/i r (log/i r a + ))(l) - (h r (\ogh r a+))(-l). 

Now notice that for b <E C°°(T) n V' + we have /i r e b = exp(h r b). With 6 = 
log/i r a + we obtain exp(/i r (log /i r a+)) = h r (h r a) = h r ia. + . Hence 

E(h r a + , h r h) = exp ((logft r 2a + )(l) - (log/i r 2a + )(-l)) . (78) 

From H Proposition 4.4] it follows that this converges to a+(l)/a+(— 1). Thus 
the proof is complete. □ 
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5 Asymptotics of the determinants of symmet- 
ric Toeplitz matrices with Fisher-Hartwig dis- 
tributions 

In this section we recall the known results about the asymptotic behavior of 
Toeplitz determinants and specialize them to the case of determinants of sym- 
metric Toeplitz matrices. Such an asymptotic formula is provided by the Fishcr- 
Hartwig conjecture, which - in the so far most general setting - has been proved 
in H (see also |l(J). We also recall the some known results about the 1Z- 
convergence of certain classes of distributions, which are later on needed in 



order to be able to use Theorem 4.2. The underlying classes of distributions are 



defined and their properties are stated next. 

For a, P, 7, S £ C and Oq £ (— 7T, it], we introduce the functions 

Lo a ,f3, go (e t9 ) = (2-2cos(6-e )) a e^ e - 9 °-*\ < 6 - 6 < 2tt. (79) 
^, eo (e i9 ) = (l-e^-^r, (80) 
^ 6o (e ie ) = (l-e i ^-^) s . (81) 

This implies 

to aM (e tf> ) = r)«+0,eo(e i0 )Z a - ,e o (e i6 ). (82) 

For a,/3,7,(5 G C with 2a £ Z_ := {-1,-2,-3,...} and 9 E (-7r,7r], 
we introduce the distributions <jJ a ,f3,e i r]-y^ and £5^0 in terms of their Fourier 
coefficients 

r , e i "('- fl °)r(l + 2a) 

K,/3,e In = 777- —3 T777- F7 — C> neZ, (83) 

I (1 + a + p — n)i (1 + a — p + n) 

JO ifn>0 
[£s,e \n - I e in(w-0 o ) ( 8 ) if n < 0. ' ^ 

It can be checked straightforwardly that if 2a ^ Z_, then 

U a ,a,0 o = 7 l2a,0 o i ^q,-q,0 o — £2a,0 o • (86) 

In what follows let GiC°°(T) stand for the set of all functions in C°°(T) 
which are nonzero on all of T and have winding number zero. In other words, 
GiC°°(T) is the set of all complex- valued functions defined on T which possess 
a logarithm that belongs to C°°(T). Moreover, let C±(T) stand for the set of 
all / 6 C°°(T) for which f n = for all n < (n > 0, resp.). We denote by 
GC£>(T) the set of all invertible functions in C±(T). 

A distribution of Fisher-Hartwig type is a distribution of the form 

C = b il W «r,/5r,0r II II ( 87 ) 

rgJWo -reAf + reM- 
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where 

(i) R > and {1, . . . , R} = Mo U M + U M_ is a decomposition into disjoint 
subsets; 

(ii) 0±, . . . , Or e (— 7r, 7r] are distinct numbers; 

(iii) b £ GiC^T); 

(iv) a r ,f3 r 6C and 2a r ^ Z_ for all r € Mo; 

(v) 7 r eC for all r e M + and <5 r e C for all r £ M_. 

The product ( j37| ) of the distributions has to be understood in the sense of (p5|). 
To this distribution we associate the function 

c(e ld ) = b(e* e ) J] u ar , 0r ,eM e ) U V^fiM 6 ) J] & r , flr (e M ). (88) 

r£Mi) reM+ reM- 

Such a function will be called a function of Fisher- Hartwig type. 
The following result has been proved in |], Proposition 5.5]. 

Proposition 5.1 Let c be the distribution ( p7| ) and c be the function (p8|). Pui 
if = {e^ : 1 < r < P}. TTien 

(a) c £ T>'(K) and c /ias i/ie smooth part c; 

(b) if M = M_ = and b £ GCf(T), iften c e gV + {K); 

(c) if M = M+ = 9 and be GC*f°(T), then c £ QV'_{K). 

Moreover, if He a > —1/2, Re 7 > — 1 and Re<5 > —1, i/ien i/ie distribution c 
can be identified with c £ L : (T) in i/ie sense of (22). 

In what follows, we agree on the following conventions. We say that 

(i) a effects ^-convergence w.r.t. [a + ,a_] and (—00, —00, —00) if and only if 
for each gel, a effects ^-convergence w.r.t. [a + , a_] and (g, g, g); 

(ii) a effects 7£-convergence w.r.t. [a+, o_] and (—00, —00, /1) if and only if for 
each gel, a effects ^-convergence w.r.t. [a + , a_] and (g, g, fi); 

(iii) a effects "^.-convergence w.r.t. [a + , a_] and (—00, /i, —00) if and only if for 
each gel, a effects "/^.-convergence w.r.t. [a + , a_] and (g, fi, g). 

Finally, let 0(N~°°) stand for a sequence of complex numbers which is 0(N B ) 
for each gel. A maximum taken oven an empty set is considered to be —00. 

Under certain conditions on the parameters, the asymptotic behavior of the 
determinants det T/v (c) with c given by ( |S7| ) is described by the Fisher-Hartwig 
conjecture. The proof of this conjecture together with the statement that such 
distributions effect 72,-convergence was the main result of (see Section 13 
therein or [|l0| Section 6]). 

Recall that Z_ := {-1, -2, -3, . . .}. Moreover, given b e GiG°°(T), let G[b] 
and E[b] stand for the constants (]?]) and (|J), and let b± e GC±(T) stand for 
the functions (Eof). 
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Theorem 5.2 Let 9i,...,6r G (— tt,tt] be distinct numbers, R > 0, put t r = 
e l9r , let{l,...,R} = M UM + UM^UM_UM* be a decomposition into disjoint 
subsets, let c be the distribution 

C = b II W «rA-,0r II Vlrfir J| & r ,0 r , (89) 

r£Mo r£M + UJK* r£M_UM* 

and assume that the following conditions are satisfied: 

(a) b G GiC 00 ^); 

(b) 2a r £ Z_, a r + (3 r <£ Z_ U {0}, a r - (3 r <£ Z_ U {0} for each r G M ; 
fcj 7 r ^ Z_ U {0} /or eac/i r G M + ; 

faV 7 r G Z_ /or eac/i r G M£; 

fej <5 r £ Z_ U {0} /or each r G M_; 

f/j 5 r G Z_ /or eac/i r G AT. 

go < (or, equivalently, Qi + Q2 < 0), w/iere 

01 = max{-l - 2Re/3 r : r G M } U { — 1 + Re<5 r : r G M_}, (90) 
£ 2 = max{-l + 2Re/? r : r G Mo) U { — 1 + Re~/ r : r G M+}, (91) 

^ 0_ \-oo if M = 9, 

g = maxjeS, £>i + p 2 }- (93) 
Finally, define the following distributions and constants: 

C+=G[b]h+ [] ^r,9r- C_=G[6]H_ n ^r,9r. ( 94 ) 

r£M UM + UM* r£M UM_UM* 

fir (95) 

r£M 

& = £wn 2(i±^S_^W jj (1 _ t , irr ,... 

rGM ° rGAfoUAf+UM* 

seM UM_UM* 

x n n mm -7 -, (96) 

r£MoUM_UM* reM UAI+UM| 

where j r = a r + (3 r and S r = a r — (3 r for r G Mo. TTien 
(i) Tjv(c) is invertible for all sufficiently large N; 
(ii) c effects TZ- convergence with respect to [c+,c_] and (go, gi, Q2); 
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(Hi) detT N {c) = G[b} N N UT E T {l + 0(N<> )). 

Since ^j^ 1 = £_i o and £_i = where £±± are the distributions defined 
in the paragraph before Proposition |4.5| , the previous theorem implies that 
the distribution h = effects ^-convergence with respect to [l,h] and 

(—00, 0, —00). This is what we have used in the proof of Theorem |4.2| . 

Now we specialize the above theorem to the case of distributions which are 
even and which have no singularities at —1 and 1. Note that 



J a r ,p r ,e r — w Qr _/3 r _e r 



*r>,0r = ilr-8r- ( 97 ) 



With this observation, it is easy to single out the class of distributions of Fisher- 
Hartwig type which are even. 

Corollary 5.3 Let 0\,... ,9 R € (-7r,0) U (0,7r) be such that \9 r \ ^ \6 S \ for all 
1 < r < s < R, R > 0, put t r = e lBr , let {1, . . . ,R} = M U M± U be a 
decomposition into three disjoint subsets, let c be the distribution 

C = b J| UJa r ,l3 r fi r 0J ar ,~l3 r -e r Yl IhrfiAlTi-BTi ( 98 ) 
reMo reM±UM$_ 

and assume that the following conditions are satisfied: 

(a) be GiC°°(T) is even; 

(b) 2a r £ a r + r <£ Z_ U {0}, a r - /3 r £ Z_ U {0} and |Re/3 r | < 1/2 for 
each r € Mq; 

(c) j r <£Z-U {0} and Rej r < 1 for each r G M±; 

(d) 7 r £ Z_ for each r 6 AQ. 
Define the numbers 

Q12 = max{-l + 2|Re/3 r : r G M } U {-1 + Re7 r : r € M±}, (99) 

*-l-±%2.tl doc) 



-00 if M = 

£»o = max{pQ, 2gi 2 }, (101) 
and the following distribution and constants: 
c + =G[b]h+ Yl 

J] »hr,9r (102) 

^ ym = 2 ^(a^-/?, 2 ), (103) 

reM 
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™ p., tt G 2 (l + a r + /3 r )G 2 (l + a r -/3 r ) 
^ = m 11 G 2 (l + 2a r ) 



x n 6 + ( ir )-%-WL(t r )- 2 t^^) n mm 



2 7 . 

re Mo r£M±L)Ml 



x n a - t r t.)" (ar_/Sp)Ca *~ /J * ) (i - 1- vr (Q - +/3 " )(Qs+/3s) 

r,s6M 



X 



J] (1- tr i-l)-(«-/Jr)<a.+/J.)( 1 _ t -l ta )-(a r +/9 r )(a.-/S.) 



r,s£M 



X 



rGM 

x n (i-^vr 7 - 73 - (104) 

r,s£ilf±UM| 

TTien 

(ij Tjv(c) is invertible for all sufficiently large N; 
(ii) c effects 1Z- convergence with respect to [c+,c+] and (qq, qyi, Qvx)', 
(in) detTjv(c) = G[fe] Ar iV n T ym £:^ m (l + 0(N«°)). 



Proof. This corollary is just a special case of Theorem 5.2. The setting in which 
we have to apply Theorem 5.2 is the following. The number R of Theorem |5.2| 
is twice the number R of this corollary. More precisely, the set Mo of Theorem 



5.2 has to be identified with two copies of the present set Mq. Let us denote 



these two copies by and M (2) . The set M+ (M*, resp.) of Theorem ^2| 

has to be identified with one copy of the present set M± (MJ, resp.). In the 
same way, the set M_ (M* , resp.) of Theorem |5.2j has also to be identified with 
one copy of the present set M± (M±, resp.). The parameters corresponding to 
the index set 

• Mq 1 ^ are {a r ,(3 r ,6 r ) for r £ Mo; 

(1) 

• Mq ' are (a r , —fir, —d r ) for r £ M ; 

• M+ U M* are (j r ,0r) for r £ M± U M|; 

• M_ U M* are (7 r , -0 r ) for r 6 M± U MjjL. 

From this it is readily seen that Qi = Q2 = Qn- In order for go, or equivalently 
£>i + Qi = 2gi2, to be negative, it is necessary and sufficient that |Re/3 r | < 1/2 
and Re 7,- < 1. H ence conditions (a)-(d) of this corollary imply conditions 



(a)-(g) of Theorem ^2 
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The distributions c + and c_ denned in (53) are given by ( |102| ) and 

c- = G[b]h- JJ U+0r,-e r U-0r,e r I] ( 105 ) 

r£M rGMj-UMJ 

Since 6 is odd, we have 6_ = 6+ and thus c_ = c+. Because of the above param- 
eters corresponding to Mq and Mq 2 \ it is easily seen that fir as given in ( p5| ) 
becomes fi^T 11 as given above. It is somewhat troublesome, but straightforward 
to verify that Et as given in (Q) becomes E^ m . This completes the proof. □ 



6 Asymptotics of determinants of symmetric 
Toeplitz plus Hankel matrices with 
Fisher-Hartwig distributions 



In this section we finally combine the Limit Theorem (Theorem 4.2) and The- 
orem 3.3 with Corollary 5.3 in order to obtain an asymptotic formula for the 



determinants of symmetric Toeplitz + Hankel matrices. This result is based on 
the asymptotic formula for determinants of (symmetric) Toeplitz matrices. 



In what follows, let G[b] and E[b] be the constants (M) and (13), and let b± 
be the functions (llC 



Theorem 6.1 Let c be a distribution that fulfills the assumptions of Corollary 
|5.3|. Define the number Q\i by (|99|) and the constants 



re M 

G{1 + a r + (3 r )G{l + a r - I3 r ) 



(106) 



EZ m = E[b] J] 



reM 



G(l + 2a r ) 



TT — - 

l k a+t 



reM (l + tr 1 )C*+W/»(l + t P )(*-«/»" (1 + ,-l)W2 

x n 6+ (t r 6_ (t r )-<° ! -+^) JJ 6_(t P )-^' 

r-GA/o reAf±UMT 

X J] (1 - t r i S )- (a ^ )K - ft)/2 (l - t -l t jl)-(«r+/Sr)(a.+/J.)/a 
r,seM 

)/2 



X 



[] (1- ir t 7 l)-(»r-/Jr)(a.+/J.)/2 (1 _ t -l ts y( ar+ r )( a3 -f3 3 



r,s£M 
r^s 



11 (l - i,T 1 <; 1 )-( a '-+ /3 "^(i - t r i s - 1 )-( Q '- /3 "^ 



rGM 
seM-tUM| 
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Then 



x n (i-v 1 *; 1 )" 7 ^ 72 - 

r,s£M±UM± 



detM N (c) = G[b] N N n M E s * m (l + 0{N S ^)) 



(107) 



(108) 



Proof. The distribution c belongs to T>'(K U K), where K — {t r : 1 < r < R}. 
Thus the assumptions of Theorem 3J3 are fulfilled and we obtain 

detT 2 Jv(crc) = {detM N (c)f. 

Moreover, the distribution c has the smooth part 

c{e w ) = b{e ie ) [] uj aTM M 9 )u ar ,-0 r ,-eM 9 ) \{ %^(e 4 %„,-0 r (e je ), 

re M r£M±UMj_ 



From Corollary |5.3| (ii) it follows that c effects ^-con verge nce w ith respect to 
[c + ,c + ] and (go, gi2, Q 12), where go is defined by (10C) and ( 101 ) and c + is the 
distribution given by (102). Notice that go and gi2 are negative real numbers. 
Obviously c + belongs to QT>' + (K U K) and has the smooth part 

c + (e* e )=Glb]h + (e M ) J] Va r+ 0Me i6 )Va r - 0r ,-e r (e ie ) R V^eM*)- 

r£M r£M±UM^ 



It is readily seen that c(t) = c + (t)c + (t). Hence the assumptions of Theorem 4.2 



are fulfilled with a and a + replaced by c and c + , respectively, and K replaced 
by K U K. Thus 

Mi) 



detT 2N (ac) =detT 27V (c) 



M-i) 



l |_ Q ^max{eo,ei2}~j 



Note that max{0o> gn} = maxj^Q, 2gi2, 012} = max{gQ, P12} = £12- Combining 
the previous formulas we obtain 



(detMAr(c)) 2 = detT 2Ar (c) 



Ml) 

M-i) 



(l + 0(N Bl2 )) 



The asymptotics of detT2jv(c) follows from Corollary 5.2 (iii). Observe the 
change from N to 2N and that go < £12 ■ Thus 



JT M-i) 



(detMjv(c)) 2 =G[b] 2N {2N) n T E s 7 



(l + 0(N ei2 )), 



where fi^ m and £y ym are the constants fllQgj ) and (|0J. Notice also that 
Ml) _ TT ,2 {a l- ^-tr r ) a ^(l-t r r^ 



~o syl 



c+(-l) 



n 22 

r£M 

MD 

X ■ 



n 



b + (-l) J-- 1 - (l + t" 1 )^ 
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Taking into account that 

Mi) 



E[bY = E[b] 



■M-i)' 

it is readily seen that 



pgr Eg* Ml) = (^ m )2. 



Obviously, f2^ m = 20^ m . Hence the last asymptotic formula becomes 

(detMjv(c)) 2 = G[b] 2N N 2(T m™ {Ef™) 2 (1 + 0{N ei2 )) . (109) 

From this the desired asymptotic formula follows, up to a sign, which will be 
determined by the following argument. Let U C C d , where d = 2|M | + \M± U 
MJ | , be the set of all d-tuples 

Z = Uar,Pr)r£M , (7r)r6A/± , ( 7 r)reA/J ) (HO) 

such that 

(a) 2a r £ Z_, a r ± r £ Z_ U {0}, |Rc/3,.| < 1/2 for all r G M ; 

(b) 7r (/ Z_ U {0}, Re 7r < 1 for all r G M± U M±. 

Let the even function 6j 6 GC°° (T) be arbitrary but fixed, and introduce 

detUjy(c) 



where c is the distribution (|9§|) with the parameters given by (110), and the 
constants and E S j^ m are defined correspondingly. 

From what has been proved so far, it follows that 

(f N (z)) 2 - (f(z)) 2 (HI) 

for z G U. Moreover, the convergence is uniform on compact subsets of U. 
Since U is connected and since f(z) / for all z G U, it follows that either 
/^f(z) — > /(z) on all of U or fjf(z) — > —f(z) on all of J7. 

Now let C7o stand for the set of all parameters (11C) for which 

(a) 2a r (f. Z_, | Re /3 r | < 1/2 for all r £ M ; 

(b) Re 7r < 1 for all r 6 M± U M£. 

Because /n(z) and /(z) are functions that depend analytically on z, because of 
the uniform convergence on compact subsets of U and because of the concrete 
structure of U and Uq, it follows that either /at(z) — > /(.z) on all of Uq or 
/ivO?) -> -/(«) on all of f/ . 

For z — we know that G f/o and that /jv(0) — > /(0). Thus we can 
conclude that Jn(z) — > /(z) on all of t/o- Hence the desired asymptotic formula 
with the correct sign follows. □ 
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7 Asymptotics of determinants of symmetric 
Toeplitz plus Hankel matrices with piecewise 
continuous functions 



From Theorem 6.1 we obtain immediately the following result concerning the 
asymptotics of the determinants of symmetric Toeplitz plus Hankel matrices 
M/v(</>) with a particular piecewise continuous generating function. 

Theorem 7.1 Let 6 G (0, n) and [3 £ C be such that |Re/3| < 1/2. Put 
t = e w ° . Then 

det M N (tg e t-p,-e a ) „ 
urn = Hi. 

where 

E = 2^(1 - tlr^d - f - 2) -4 (l-*o-^(l-*o)-; G(1 + mi _ p) , 

(l + *o ) 2 ( 1 + i o) 2 



Proof. In Theorem |6J we put M± = MJ = 0, M = {1}, a.\ = 0, ft = [3, 
6i = 6 and b(t) = 1. Observe that wo,/3,0 o w o,-/?,-0 o = t/3,e t-p,-e - D 

The main result concerning the asymptotics of the determinants of sym- 
metric Toeplitz 4- Hankel matrices M/v(</>) with "general" piecewise continuous 
generating functions is the following theorem. 

Theorem 7.2 Let 

R 

c{e w ) = b(e* e ) Y[t^e r (e w )t-P r ,-e r (e l9 ) (112) 

r=l 

where b G G\B\ is an even function, 6\,...,6r G (0, tt) are distinct numbers, 
and /3t, . . . , (3r £ C are suc/i that \Ke f3 r \ < 1/2 /or all 1 < r < R. Let G[b] 
and E[b] be the constants (Jfy and (13), b± be the functions (0), t r — e l9r , 
1 < r < R, and introduce the constants 

" S // m = -X> 2 , (113) 

r=l 

£7 = J] G(l + &)G(1 - (3 r )(l - tl)-^ 2 (l - t; 2 )-^/ 2 

r=l 

x n a - t r t s )-^°(i - t- vr^a - M7 x )^«(i - t-H s f^ 

l<r<s<R 

x TT 2-A 2 ii^-i_ii^ TT 6 + (t P )^6_(t r )-^. (114) 
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Then 



detMjv(c) 



(115) 



Proof. The asymptotic formula follows from Corollary 2.5), Theorem 7.1 and 
the asymptotic formula (|l2|). □ 



We remark that in the case b € GiC°°(T), the previous theorem follows also 



directly from Theorem S.l 



In view of the asymptotic formulas established in Theorem |2.3| , Theorem p.4| 
and Theorem 7.1 we now raise the following conjecture. 

Conjecture 7.3 Let 6q £ (0, n) and Pi,ffa & C be such that Re/3i| < 1/2 and 
|Re/3 2 | < 1/2 and \Re f3 1 +Re/3 2 | < 1/2. Put t = e w « . Then 



where 



,. dct M N (tp lt g ti3 2 -e a ) _ 



(116) 



£7 = G(l + /3i)G(l + /3 2 )G(1 - A - /3 2 )2 ft& 

X(l - t- 2)/5 1 2 /2+ft& (1 _ f 2 ) /3 2 2 /2 + ^ 1 ^ 2 

(l-t-yi/a^-tojA/a 



(117) 



(118) 



(l+to^^a+to)*/ 2 ' 

This conjecture fits with the results established in Theorem 2.S, Theorem |2.4| 



and Theorem 7.1 
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